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Abstract: We compute all two- loop master integrals which are required for the evaluation of next-to- 
leading order QCD corrections in Higgs boson production via gluon fusion. Many two-loop amplitudes 
for 2^1 processes in the Standard Model and beyond can be expressed in terms of these integrals 
using automated reduction techniques. These integrals also form a subset of the master integrals for more 
complicated 2 — > 2 amplitudes with massive propagators in the loops. As a first application, we evaluate 
the two-loop amplitude for Higgs boson production in gluon fusion via a massive quark. Our result is the 
first independent check of the calculation of Spira, Djouadi, Graudenz and Zerwas. We also present for 
the first time the two-loop amplitude for gg ^ h via a massive squark. 



1. Introduction 



Next-to-leading (NLO) QCD corrections are important for a wide range of processes at the LHC. Recently, 
there have been new breakthroughs in developing efficient techniques for one-loop amplitudes, e.g. [1]. It 
is realistic that these theoretical advances will improve the data analysis for many interesting observables. 
However, processes which cannot occur with tree-level interactions at leading order, require two-loop rather 
than one-loop amplitudes at NLO. 

Loop induced processes should be rather sensitive to new physics (see for example [2]). It is possible 
that we revisit their NLO corrections several times in the future, in order to include new types of particle 
interactions in the loops. In this paper, we study a basic LHC process of this kind: the production of a 
Higgs boson in gluon fusion. 

The two-loop QCD corrections for gg ^ H in the minimal Standard Model and its two-Higgs-doublet 
extensions have already been computed in [3] . In this calculation, the mass effects of the quark coupled to 
the Higgs boson were also fully accounted for. The impact of the NLO correction is striking; for example, 
the SM Higgs boson production cross-section increases by more than 70%. 

This calculation has never been verified in the literature, except in well known limits such as within 
the heavy top-quark approximation [4]. It is thus important to perform an independent computation. 
Our main objective, however, is to automatize the evaluation of the two-loop amplitude in QCD for this 
and other processes with similar distribution of massive particles in the loops. This is essential for future 
applications in gluon fusion and other processes for a variety of extensions of the Standard Model. 

In 2001 Laporta introduced a new algorithm for an automated reduction of multi-loop integrals to 
master integrals [5]. A parallel rapid development of the Mellin-Barnes method [6-10], the differential 
equation method [11], and sector decomposition [12-14], yielded robust technology for computing master 
integrals. As a result, two-loop calculations for three and four point functions with more than one scale 
are now tractable. 

In this paper we study the two-loop integrals which are required in gluon fusion processes. We apply 
the algorithm of Laporta, using the package AIR [15] and an independent MATHEMATICA implemen- 
tation [16], to perform a reduction to master integrals. Some master integrals were already known in 
the literature [17-23]. We have recomputed them using the method of differential equations [11]. We 
present here for the first time the remaining master integrals, including the most complicated scalar cross- 
triangle. The same master integrals enter the evaluation of two-loop amplitudes of more complicated 2^2 
processes, such as heavy quark pair production. 

Our results are first given as an expansion in the dimension parameter e = (4 — (i)/2 in terms of har- 
monic polylogarithms. We present the series coefficients through the order where harmonic polylogarithms 
with transcendentality four appear. The polylogarithms are real valued in a non-physical kinematic region 
corresponding to imaginary center of mass energy s < 0. Then we perform explicitly the analytic continua- 
tion in the physical regions below and above the threshold corresponding to two on-shell heavy particles in 
the loops: s = 4mf. For s > 4mf the analytic continuation proceeds as in [24,29]. For s < 4m^, we find the 
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analytic continuation of harmonic polylogarithms using the procedure described in [19]. Interestingly, the 
result for the master integrals is in general simpler than the analytic continuation for individual harmonic 
polylogarithms. 

As a first application, we compute the two-loop amplitude for gg ^ h in the Standard Model. We 
have compared our result with the expression in [25], in the non-physical region. The result of [25] was 
derived by series expanding the integral representation of the two-loop amplitude from [3] in a kinematic 
variable, and mapping the expansion to a carefully chosen ansatz. Our result agrees with [25]. This is the 
first independent check of the two-loop amplitude in [3]. 

We present here a new result for the two-loop amplitude gg ^ h via scalar-quarks. In the heavy 
squark limit, our result agrees with [27]. With a completed setup for the reduction procedure and the 
expressions for the master integrals, the evaluation of this new result is fully automated. Other amplitudes 
with different particle content in the loops can be obtained easily. We note that preliminary numerical 
results for the NLO K-factor for the squark case have been presented in [28] . 

2. Reduction of the amplitudes 

We consider QCD virtual corrections to the process 



with just one massive particle running in the loops. Some typical Feynman diagrams are shown in Figure 
1 for the cases of a heavy quark and a squark in the loop. 

It is convenient to project the two-loop amplitudes onto scalar form factors. In this way we are left 
only with loop integrals involving scalar products in the numerator. The scalar integrals can be classified 
into topologies according to their denominators. In gg ^ h we find diagrams with at most six propagators, 
of which five can be massive. In the numerators, however, one can find seven independent scalar products. 
The irreducible scalar product can be dealt with by introducing an additional propagator, which is raised 
to negative powers in the expressions for the physical amplitudes. After the introduction of the auxiliary 
propagators, all scalar integrals belong to subtopologies of the three topologies, shown in Figure 2, with 
denominators: 



g + g ^ h 
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Figure 1: Typical Feynman diagrams in the two loop contributions to gg H with (a) a heavy fermion in the 
loop, (b) a heavy scalar in the loop . 



where pi2 = P1+P2, "T-t is the mass of the particle running in the loops and k and I are the loop momenta. 

The reduction to master integrals is done using integration by part identities [30, 31] combined with 
the Laporta algorithm [5] in [15, 16]. We found 17 master integrals, which are shown in Figure 3. It is 
possible to choose a different basis of master integrals; the basis we choose is particularly convenient for 
the method of differential equations. 

The master integrals in the first two lines of Figure 3 are products of known one-loop integrals [17, 19]. 
The master integrals in the third, fourth and fifth line in Figure 3 are non-factorizable. Integrals in the 
third and fourth line were calculated already in [18]^ and [19,21,22]. respectively. The double triangle, 
last diagram in the third line was calculated in [21-23]. Also the six propagators triangle - third diagram 
in the last line of Figure 3 - has been calculated in [20]. 

3. Master integrals 

We computed all master integrals using the differential equation method [11,32-36]. The natural variable 
to express the results is 

VT^^ - 1 . , 4mt^ 



+ ie where r = , (3-1) 



^Our results fully agree with the results quoted in this reference taken from the electronic file in 
http://pheno.physik.uni-freiburg.de/bonciani/. The printed version contains several typographical mistakes. 
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Figure 2: Master topologies. Wavy lines denote massless particles both external and internal. Internal massive 
lines arc denoted by single straight lines whereas the double line denotes a massive external particle. 

with s = {pi +P2)'^- The variable x is real valued in the space-like region (s < 0) and in the physical region 
above threshold {s > Ami). Below threshold x lies on the unit circle in the complex plane. In that region, 
we introduce the variable 9 such that x = e^^ . For quick reference, in Table 1 we list the domain of each 
variable in the different kinematical regions. 



Region 


s 


T 


x 


space- like 


-00 < s < 


> T > -00 


< X <1 


below threshold 


< s < Amt^ 


00 > T > 1 


X = e*^ with < 6* < vr 


above threshold 


Arrit^ < s < 00 


1 > r > 


-1 < X < 



Table 1: Domain spanned by the variables in the different kinematical regions 

The dependence on x of the master integrals is determined by solving the associated differential equa- 
tions. These are obtained by taking the derivative with respect to x of the loop integrals and exchanging 
the order of the differential operator and the loop integration. In this way, the derivative of a given master 
integral is expressed in terms of scalar integrals which can again be reduced to masters. Applying this 
procedure to all integrals in our master basis, we derive a closed system of differential equations. 
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Figure 3: Set of master integrals. The conventions for the lines are as in Figure 2. Each dot on a propagator line 
denotes an additional power of the propagator in the denominator. The diagram with a big dot contains a numerator, 
it is defined in the following section. 

We solve the differential equations order by order in powers of e. Only integrations with kernels l/x, 
1/(1 — x) and 1/(1 + x) are required, and the solutions can be written exclusively in terms of harmonic 
polylogarithms [29]. 

To fully determine the solution of the differential equations, we require the value of the master integrals 
at a certain kinematic point. The value at x = 1 is very easy to obtain. This limit corresponds to setting 
the external momenta to zero (s = 0). All non-factorizable master integrals MI^^^^ collapse to a vacuum 
sunset diagram with extra powers of propagators: 



lim Ml(NF) = 




1^3 



The exponent f2 corresponds to the number of massless propagators in the integral whereas vi + 1^3 is the 
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number of massive propagators. One can easily compute: 



1^1 

We have observed that one could fix the solution of the differential equations by requiring simply that the 
X — > 1 limit is finite, since the homogeneous solutions usually diverge at x = 1. The explicit formula for 
the limit x = 1 in Eq. 3.3 was then an additional consistency check of our calculation. 

In one master integral, the x = 1 limit does not commute with the expansion around e = 0, due to a 
coUinear singularity as s vanishes. For this master integral, we have used the massless limit x — > 0, which 
is well behaved: 



with 




In the following sections we shall present our results in the space-like and < s < 4m^ regions, and de- 
scribe the analytic continuation procedure. We have performed several checks on our results for the master 
integrals. We have verified that our expressions satisfy the differential equations before and after analytic 
continuation. We have also computed all master integrals numerically from their Feynman parameteriza- 
tion, by using sector decomposition [12-14]. Our analytic expressions agree fully with the direct numerical 
evaluation. 



4. Results in the region s < 

We now present the results for the master integrals in the space-like region. We give the definition of the 
master integrals in terms of the propagators listed in section 2, and their e-expansion in terms of harmonic 
polylogarithms. 

4.1 One loop integrals 
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7 ^TT^/^ D21IJ23 1-e 



^mbub(^) = 1 (4.4) 

i^mbub(^) = ^{-x + (x + l)F(0;x) + l} (4.5) 
Fibub(2;) = ^{^(-^'x-12x-7r2 + 12) +(x + l)F(0;x)-2(x + l)i7(-l,0;x) 

+ (x + l)F(0,0;x)} (4.6) 
i^^bub(^) = :r^{-^vr2(x + l)-2((2 + C(3))x + C(3)-2) 

1 — X 

+ Itt\x + l)Hi-l;x) - I (-12 + TT^) {x + l)i7(0; x) 

- 2{x + 1)H{0, -1, 0; x) - 2(x + 1)H{-1, 0; x) + (x + 1)H{0, 0; x) 

+ 4(x + 1)H{-1, -1,0; x) - 2(x + 1)F(-1, 0, 0; x) + (x + 1)H{0, 0, 0; x)| (4.7) 

i^^bub(^) = ^^'(^ + 1) - l^'i^ + 1) - 2((4 + C(3))x + C(3) - 4) 

+ ^tt\x + 1)H{0, -1; x) + ^(x + 1)H{-1; x) (vr^ + 12C(3)) 

- ^(x + 1) (tt^ + 12(-2 + C(3))) H{0; x) - 2(x + 1)H{0, 0, -1,0; x) 

- 2(x + 1)F(0, -1, 0; x) - '^tt\x + l)i^(-l, -1; x) 

+ I (-12 + vr^) (x + 1)H{-1, 0;x)-l (-12 + vr^) (x + 1)H{0, 0; x) 
3 b 

+ 4(x + 1)F(0, -1,-1,0; x) -2(x + 1)F(0, -1,0,0; x) + 4(x + 1)F(-1, 0, -1, 0; x) 
+ 4(x + l)i^(-l, -1,0; x) - 2(x + 1)F(-1, 0, 0; x) + (x + 1)H{0, 0, 0; x) 
+ 4(x + 1)F(-1, -1,0,0; x) -8(x + l)i7(-l, -1,-1,0; x) 

+ (x + 1)H{0, 0, 0, 0; x) - 2(x + 1)-H'(-1, 0, 0, 0; x)| (4.8) 
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6 

+ 2F(0, -1, 0, 0; x) + F(0, 0, 0; x) - i7(0, 0, 0, 0; x)} (4.12) 
4.2 Factorizable integrals 
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4.3 Three propagator integrals 
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+ 24H{0, -1, 1, 0; x) + 12H{0, 0, 0; x) - 2AH{1, 0, -1,0; x) - 8H{1, 0, 0; x) 

+ 8H{1, 0, 1, 0; x) + 24H{0, 1, -1, 0; x) - 20H{0, 1, 0, 0; x) - 8H{0, 1, 1,0; x) 

- UH{0, 0, 0, 0; x) + 12F(1, 0, 0, 0; x) - 8H{1, 1,0, 0; x)| (4.23) 
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+ ^(x - 1)H{1; x) (6C(3)(7 - 4x) + tt\x - 1)) 

- \tt^{x - l){5x - 3)H{0, 1; x) - 6{x - l){13x - 7)H{0, 0, -1, 0; x) 

+ 12(x - l)(5x - 3)H{0, -1, 0; x) - 67r2(x - 1)'^H{-1, -1; x) 
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2 

- -tt'^H{1, 1; x)(x - 1)2 - 4(5x - 3)if(0, 1, 0; x){x ~ 1) 

+ 2(x - l)(13x - 7)H{0, 0, 1, 0; x) + 36(x - l)(5x - 3)H{0, -1, -1, 0; x) 

- 24(x - l)(5x - 3)H{0, -1, 0, 0; x) - 12(x - l)(5x - 3)H{0, -1, 1, 0; x) 
+ 144(x - l)^i?(-l, 0, -1, 0; x) - 72(x - lfH{-l, -1, 0; x) 

+ 48i?(-l, 0, 0; x)(x - if + 24F(-1, 1, 0; x)(x - 1)^ 

- 48i?(-l, 0, 1, 0; x)(x - if - 2(13x - 7)H{0, 0, 0; x)(x - 1) 
+ 24(x - lfH{l, -1, 0; x) - 12(x - l)(3x - 5)i7(l, 0, -1, 0; x) 

- 8i7(l, 1, 0; x)(x - if - 4(3x - 5)iJ(l, 0, 0; x)(x - 1) 

+ 4(x - l)(3x - 5)H{1, 0, 1, 0; x) - 12(x - l)(5x - 3)H{0, 1, -1, 0; x) 
+ 2(x - l)(27x - 17)if(0, 1, 0, 0; x) + 4(x - l)(5x - 3)H{0, 1, 1, 0; x) 
+ 144(x - lfH{-l, -1, 0, 0; x) - 216(x - lfH{-l, -1, -1, 0; x) 
+ 72(x - lfH{-l, -1, 1, 0; x) - 60(x - lfH{-l, 0, 0, 0; x) 
+ 72(x - lfH{-l, 1, -1, 0; x) - 48(x - lfH{-l, 1, 0, 0; x) 
+ (x - l)(29x - 15)if(0, 0, 0, 0; x) - 24(x - lfH{-l, 1, 1, 0; x) 
+ 72(x - lfH{l, -1,-1, 0; x) - 48(x - lfH{l, -1, 0, 0; x) 
+ 2(x - l)(7x - 13)H{1, 0, 0, 0; x) - 24(x - lfH{l, -1, 1, 0; x) 
+ 4(x - l)(5x - 3)i7(l, 1, 0, 0; x) - 24(x - 1)^-H"(1, 1, -1, 0; x) 

+ 8(x- l)2i7(l, 1,1,0; x)} (4.28) 
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4.4 Four propagator integrals 




d'^k f 1 



J DuDuDi^Dn 

2 1 

(^Y^) {m',)-'^Ee%ix) + 0{e') (4.29) 



F,-\x) = - (4.30) 



F,-\x) = -\ (4.31) 
F^{x) = .^^-i-^|-3x2 + (6-4C(3))x + 2(x2-l)F(0;x)-3 

- F(0, 0; x)(x - if + 2xi7(0, 0, 0; x) + 4xH{l, 0, 0; x)| (4.32) 

^3 (^) = + C(3))x2 - 2(-12 + C(3))x + ^ - -vr^ {x' - l) + 3(-4 + C(3)) 

+ ^H{0; x) {TT^ix - 1)2 + 12 (4x2 - 3C(3)x - 4)) - 12xH{l; x)C(3) 
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F^\x) = \ (4.35) 
F^\x) = ^^^{^(x-l)2 + (l-x2)F(0;x)} (4.36) 

f2{x) = -^^-^{5x2 + 2(-5 + 2C(3))x+ (^-3x^ + ^ + 3^ i/(0;x) + 5 

+ 2 (x^ - l) i?(-l, 0; x) - (x - l)(x + 2)F(0, 0; x) + (l - x^) F(l, 0; x) 

+ 2xi7(0, 1, 0; x) + xif(0, 0, 0; x)| (4.37) 

Fl{x) = ^^^{-(-16 + C(3))x2+(-32-H^ + C(3))x-4(-4 + C(3)) 

+ i?(0; x) (^-lOx^ - ^7r2(x + l)x - 3C(3)x + 10^ - '^'k'^xH{{), -1; x) 

+ ^7r2x//(0, 1; x) + \h{1- x) (-48C(3)x - vr^ (x^ - l)) 

- 2xiJ(0, 0, -1, 0; x) + 2(x - l)(x + 2)F(0, -1, 0; x) + 6 (x^ - l) H{-1, 0; x) 
+ (^-3x2 + i (-18 + tt'^)x + 6^ H{<d, 0; x) + (^-3x2 - + 3^ i/(l, 0; x) 

- 2 (x^ + X - 1) i?(0, 1, 0; x) + 2xi?(0, 0, 1, 0; x) - 2xiJ(0, -1,0, 0; x) 

- 4xi7(0, -1, 1, 0; x) + (4 - 4x2) H{-1, -1, 0; x) + 3 (x^ - l) i?(-l, 0, 0; x) 
+ 2 (x^ - l) i?(-l, 1,0; x) + (-2x2 - 3x + 4) F(0, 0, 0; x) 

+ 2(x2-l)/7(l,-l,0; x) + (— x^ — 2x + 3) H{1, 0, 0; x) 

+ (2 - 2x2) F(l, 1,0; x) - 4xi7(l, 0, 1, 0; x) - 4xi7(0, 1, -1, 0; x) 

+ 4xiJ(0, 1, 0, 0; x) + 4xiJ(0, 1, 1,0; x) + 3xi7(0, 0, 0, 0; x) 

- 2xiJ(l, 0,0,0; x)| (4.38) 
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F^(x) = 3__|;^(31x^-140x^ + (218-64C(3))x2-140a; + 3l) 

(1 — xj^x L 32 ^ 



+ ^ (-15x^ + 54x^ - 47r2x2 - 54x + 15) i7(0; x) 
+ J (x^ - 4x^ + 4a; - l) 0; x) 



1 / 1\ 

+ - {-x"^ + ix^ + 3^2 - 8x + 2) i?(0, 0; x) + ( -— + x^ -x + -j H{1, 0; x) 

- H{0, 1, 0; - ^H{0, 0, 0; x)^^ | (4.42) 

FUx) - 1 f 1^189 C(3)\ 4 I ( 233 tt^ A 3 /743 Uvr^ 3C(3)\ , 

^^^"'^ - (l-x)2xlU4 4r+lv 16 + 24 + ^^^V V 32 + 180 4 J 



x) 



1 189 1 

- - (699 + 27r2 - 192C(3)) x - C(3) + + 37r'^'i^(0, -1; x) 

+ -^^(0; x) (- (93 + 27r2) x^ + (342 + Sn^) x^ + 6 (tt^ + 12C(3)) x^ - 342x + 93) 

+ ^H{l;x) (96x^^(3) - tt^ (x^ - 4x^ + 4x - l)) 

+ x^H{0, 0, -1, 0; x) - ^Tr'^x^H{0, 1; x) 

+ ^ (x^ - 4x^ - 3x2 + 8x - 2) iJ(0, -1, 0; x) 

+ ^ (Sx^^ - 18x2 + 18x - 5) H{-1, 0; x) 

+ f-^ + 2x=^ - i- (-57 + 27r2) x^ - 5x + /f(0, 0; 
y 8 24 4 J 

+ ^ (-15x^ + 60x2 ^ g^2^2 _ ^ -^5^ ^^^^ Q. 

+ i (-x"^ + 4x2 + 3x2 _ 43. + 1) ^(0, 1, 0; x) - x^H{0, 0, 1, 0; x) 

+ iJ(0, -1, 0, 0; x)x2 + 2ii"(0, -1, 1, 0; x)x2 + (-x^ + 4x2 - 4x + l) H{-1, -1, 0; x) 

+ ^ (x^ - 4x2 + 4x - 1) H{-1, 0, 0; x) 

+ i (.T^ - 4x2 + 4x - 1) H{-1, 1, 0; x) 

/ x"' 9x2 \ 

+ f-y + 2x2 + — -4x + ljff(0,0,0;x) 

+ ^ (x^ - 4x2 + 4x - 1) H{1, -1, 0; x) 
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+ 2H{1, 0, 1, 0; x)x^ + ^ (-x^ + - 4x + l) H{1, 1, 0; x) 
+ 2H{0, 1, -1, 0; - 2H{0, 1,0, 0; x)x2 - 2H{0, 1, 1, 0; x)x'^ 

+ x2iJ(l, 0, 0, 0; x) - ^x2if(0, 0, 0, 0; x)| (4.43) 



F^(x) 



F^(x) 




J DuDuDiaDf^ 
^ i=0 



xH{0,0;x) 



9C(3) 



2(x- 1)2 

: (1 ^^)2 { - ^^'-^-(0; x) - H{0, -1,0; x) - H{0, 0; x) ^ 

- /^(0, 1, 0; x) + ^H{0, 0, 0; x) + H{1, 0, 0; x)} 

X r 1 1 llvr^ 

: -vr2i?(0, -1; x) + -vr2i/(0; x) + 9C(3) + ^ 

+ 3C(3)ff(l; x) - -7r^H(0, 1; x) - H(0, 0, -1,0; x) + 2H(0, -1, 0; x) 
6 

+ ^ (6 - vr^) H{0, 0; x) + ^vr2i?(l, 0; x) + 2H{0, 1,0; x) 

+ 2H{0, -1, -1, 0; x) + 2H{0, -1, 1, 0; x) - H{0, 0, 0; x) - 2H{1, 0, -1,0; x) 

- 2H{1, 0, 0; x) + 4i7(l, 0, 1, 0; x) + 2H{0, 1, -1, 0; x) - H{0, 1, 0, 0; x) 

- 2H{0, 1, 1, 0; x) + ^H{0, 0, 0, 0; x) + 4H{1, 0, 0, 0; x) + 2H{1, 1, 0, 0; x)| 
4.5 Five propagator integrals 



(4.44) 
(4.45) 



(4.46) 



(4.47) 




r d^k r 
Jj^ J' 



d/2 

r(i + e) 

1 -e 



ITT"-/^ J ?7r'^/2 D22D23D24D26D27 
2 

2N-2e-l 77O/ 



imlr'^-'F?ix) + Oie') 



(4.48) 
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4 

TT X 



(1-X)2 

- a;i?(0, 0, 1, 0; x) - 2xH{l, 0, 1, 0; x) - 2xH{0, 1, 0, 0; x) 
-4xiJ(l,l,0,0;x)} 



3xH{l,0,0,0;x) 



(4.49) 



Fiix) 




J' 

T(l + e 



1 



?^'^/2 DuDisDuDiaDn 

2 1 



^ ^ i=0 



I - 2i?(0, 0, 1; x) - 2H{0, 1, 0; x) + 4H{1, 0, 0; x) - 6C(3)| 
12C(3)F(0; x) + ^7r2i/(0, 1; x) - 24H{1- x)C(3) - ^ 



8F(0, 0, 0, 1; x) - 10i?(0, 0, -1, 0; x) + 4F(0, -1, 0, 1; x) - -7r'^H{l, 0; x) 



(4.50) 



(4.51) 



- 4i?(l, 0, 0, 1; x) - 4F(0, 1,0, 1; x) - 4i7(0, 0, 1, 0; x) - 4i?(0, 0, 1, 1; x) + 4ii"(0, -1,0, 0; x) 
+ 4i7(0, -1, 1,0; x) - 24F(1, 0, -1,0; x) + 4F(1, 0, 1, 0; x) + 4F(0, 1, -1,0; x) 

- 6F(0, 1,0, 0; x) - 4F(0, 1, 1, 0; x) + 12F(1, 0, 0, 0; x)| (4.52) 
4.6 Six propagator integrals 




d'^k f d'^l 1 

2 



d/2 / 

r(i + e 

1 - e 



2\--2e-2 Z70 
t . 



F^{x) + 0{e'] 



(4.53) 



(1 - X)3(x + 1) 



|8C(3)F(0; x) + 16if(0, 0, -1, 0; x) + 



vr 
10 



+ -7r2i/(0, 0; x) - 4H{0, 0, 1,0; x) - 8H{0, -1,0, 0; x) + 14F(0, 1,0, 0; x) 
3 



+ 



F(0,0,0,0;x)} 



(4.54) 



- 16 - 




2 



i= — l 

FroHx) = ^^^^{-^7r2ii'(0;x)-8F(0,-l,0;x)+4F(0,0,0;x)-12C(3)} (4.56) 

Ki-) = (r^{r'^(0'-i'-) + 24c(3)-^ 

+ ^ (vr' - 33C(3)) H{0; x) - ^n^H{{), 1; x) - 48F(1; x)C(3) 

10 8 

- 56i7(0, 0, -1, 0; x) + 167^(0, -1,0; x) - —-K^HiO, 0; x) - -7t'^H{1, 0; x) 

o o 

+ 8F(0, 0, 1, 0; x) + 64i?(0, -1, -1,0; x) - 40H{0, -1,0, 0; x) - 16i?(0, -1, 1, 0; x) 

- 8H{0, 0, 0; x) - 32H{1, 0, -1, 0; x) - 16H{0, 1, -1, 0; x) + 8F(0, 1, 0, 0; x) 

+ 12H{0, 0, 0, 0; x) + 16F(1, 0, 0, 0; x)} (4.57) 

5. Analytic continuation 

The expressions in section 4 correspond to the unphysical, space-like region. The results must be 
analytically continued towards the time-like region. Due to the threshold in s = 4mj, the physical region 
splits into two subregions, namely above and below threshold. 

5.1 Analytic continuation above threshold 

The region above threshold corresponds to the range — 1 < x < 0. The analytic continuation to 
this region is straightforward, using the properties of harmonic poly logarithms under the transformation 
X — > —X [24, 29]. For harmonic polylogarithms H(an, . . . , ai; x + ie) with ai ^ the analytic continuation 
is obtained trivially 

H{an, ...,±l;x + ie) = (-l)±i+-+'^"if(-a„, . . . , ^1; -x) . (5.1) 

where = —1,0,1 for k ^ 1. We can eliminate higher rank harmonic polylogarithms with ai = by 
applying integration by parts and product identities [29]. For instance: 

H{1, 0; x) = H{0; x) H{1; x) - H{0, 1; x) . (5.2) 
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Recursively, one can write similar identities for the harmonic poly logarithms of higher rank. At the end of 
this procedure we only require the analytic continuation of simple logarithms: 



H{l;x + is) = -H{-l;-x + i£) 
H{0; x + ie) = H{0; -x + ie) + in 
H{-l;x + ie) = -H{l;-x + ie) . 



(5.3) 



The analytic continuation described here is incorporated in the Mathematica package [37] . The expressions 
for the master integrals above threshold can be easily obtained from the ones in the space-like region using 
the routines implemented in this package. 

5.2 Analytic continuation below threshold 

Below threshold, the variable x lies on the unit circle of the complex plane. For this analytic continu- 
ation we follow the procedure in [19], which we summarize here. 

We first express our results in terms of the variable 9 given by x = exp(i 6), and introduce the following 
notation for the harmonic poly logarithms as functions of 6: 



Hc{an, ...,ai;0)'^= H{an, ai; e*^) 



(5.4) 



We now eliminate polylogarithms Hc{an, . . . ,ai;6) with a„ = 1, using integration by parts and product 
identities [24]. Then we use the analytic continuation of harmonic polylogarithms of weight one as kernels. 
For < < TT we have 



Hc{l;9) = -ln2 
i?c(0; 9) = i9 
H,{-l-9) = In 2 



sm ■ 



cos- 



+ i 



+ ^2 



vr 9 



(5.5) 
(5.6) 
(5.7) 



We can find the analytic continuation for the harmonic polylogarithms of higher weights recursively, using 



Hc{an, a„_i, ...,ai;9) = H{an, an-i, ...,ai-l) + i / d9'g{an\ 9')Hc{an-i, • • • , ai-9') 

Jo 



where 



9(0; 9) 



1 - e*' 



1 1 9 

-— + i— cot — 

2 2 2 



ei9 1 ' 

e'^ 1 19 



(5.8) 

(5.9) 
(5.10) 
(5.11) 
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In this way, we have obtained expressions for the analytically continued harmonic polylogarithms through 
weight 4. The analytically continued master integrals can be expressed in terms of the following functions 
[19,38,39]: 



ch{e) 



Lscij(^) 



In 2 



sm — 
2 



sm ■ 



re 

/ d^'ln^-^: 
Jo 



sm ■ 



2 

In-'-i 2 



d' 



cos ■ 





. 9' 




9' 




sin — - 


2 


cos — 




2 




2 



(5.12) 
(5.13) 
(5.14) 
(5.15) 



Interestingly, these functions are a smaller set than the functions that appear in the analytic continuation 
of individual harmonic polylogarithms. 



6. Results for the master integrals in the region below threshold 

We now present the master integrals in the kinematic region < s < 4m|. The expressions Fj{9) 
correspond to the analytic continuation of the coefficients Fj(x) in section 4. 

6.1 One loop integrals 



^-1 

mbub 



(9) 



^mbub(^) = 



mbub 



(9) 



^mbub(^) 



-^COt I 2 I +1 



■- cot ( - ^ - 29Ch{9 - tt) + 2Ls^°V - tt)) + 2 
^mbub(^) = I cot (Jj (6( - 20 - ^) - 129Ch{9 - tt) - 129Ch{9 - 
+ (24Cli(e - tt) + 12)Lsf V - tt) + 12Ls^°^(0 - tt)) + 4 



^ cot ( - TT^ - 240 + 6^ - 40 - y)cii(0 - vr) - 120Cli(0 - nf 

- 80Cli(0 - nf + (24Cli(0 - + 24Cli(0 - tt) + 24) Lsf\9 - tt) 
+ (24Cli(0 - tt) + 12)Ls^°^(0 - tt) + 8Lsi°^(0 - tt) + 127rC(3)) + 8 



(6.1) 
(6.2) 

(6.3) 



(6.4) 



(6.5) 
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fLAO) = (6-6) 
= I csc^ (e' + {AO - 87r)Ls(°V - ^) - SLs^^V - ^) - 14C(3)) (6.7) 

+ 12^Ls^°^ - tt) + 24Ls^^^ - tt) + 42C(3)) (6.8) 
6.2 Three propagator integrals 

m = (6-9) 

fI (9) = ^ csc^ (^^^ (d^ + ^^Cli (9) + 2eLsf (e) + 6Lsf (^ - tt) (-27r + ^) 

- 6Ls|^^(^) - 12Ls^^^(0 - tt) - 21C(3)) (6.10) 
F^{9) = csc^ ( ^ ) hdn* - ISOvr^^ - 60^^ - 30^^ - l209^Ch{9f + UOtt"^ log2(2) - 1201og^(2) 

- m9Lsf\9) + A80Lsf\9f + Lsf\9 - 7r)(28807r - 14400 + 2880Ls^°^ (0)) 

+ 216OLs5°^(0 - Trf + 2A09Lsf\9) + 720Lsf ^(0 - tt) ( - tt + 30) + 144OLs^^^(0) 

+ 2880Ls^^V - tt) + 18OLs?^(20) - 720LsJV - tt) + 144O0Lsc2,2(0) + 144OLsLsci,i,2(0) 

- 2880Li4 + 5040C(3) - 25201og(2)C(3) + Cli(0)( - 2400^ - 4809LsP (9) 

- 1440(0 - 27r)Ls^°^(0 - tt) + 144OLs^^^(0) + 288OLs|^^(0 - tt) + 5040C(3)) ) (6.11) 



^2"'(^) = ^^csc(0) (6.12) 



F2°(0) = (.Q^ ^ j gsg(^^ + ps^.(^^) - + 0Cli(0) + 30Cli(0 - tt) - Ls^°^(0) - 3Ls^°V - ^r)) 

(6.13) 

(0) = + ^ cot Q csc(0) ( - 02 - 02cii (0) - 20Ls^°^ (0) + (127r - 60)Ls^°^ (0 - tt) + 6Ls^^^ (0) 

+ 12Ls^^^(0 - tt) + 21C(3)) +^^^(^ + 37r^ + 20 + 20^ + 40Cli(0)2 - 240Cli(0 - tt) 
+ 360Cli(0 - Trf + Cli(0)(240Cli(0 - tt) - 80) + (-8Cli(0) - 24Cli(0 - tt) + 8)Ls^°^(0) 
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+ (-24Cli(0) - 72Cli(0 - tt) + 24)Ls^°V - ^r) - 4Ls^°V) - 36Lsf ^(0 - tt) 

- 24Lsc2,2(0)) (6.14) 
F^{e) = cot (-) csc(0)(- 797r^ + 18O7r30 + 60^2 ^300^ + 120^2^1^(^)2 

- 1207r2 log2(2) + 1201og^(2) - 480Ls^°V)^ + l^sf\e){A80Ch{9)e + 4806 

- 2880Ls^°V - ^)) + (240(60 - 127r) + 1440(0 - 2^)Cli(0))Ls^°V - tt) - 2160Ls5°V - nf 

- 2409Lsf\e) + 240(3^ - 90)Ls^°^(0 - vr) + (-144OCli(0) - 1440)Ls^^^ (0) 

+ (-288OCli(0) - 288O)Ls^^^(0 - tt) - 18OLsJ^(20) + 72OLsJ^(0 - tt) - 144O0Lsc2,2(0) 

- 144OLsLsci,i,2(0) + 2880Li4 + 24OCli(0) (0^ - 21C(3)) 

- 5040C(3) + 25201og(2)C(3)) + (40Cli(0)^ + (180^ + 180 + 277r3) Cli(0 - tt) 

- lO80Cli(0 - irf + lO80Cli(0 - vr)^ + Cli(0)2(360Cli(0 - tt) - 120) 

+ Cli(0) (603 + lO8Cli(0 - Trfe - 72Cli(0 - 7r)0 + 60 + Qtt^) + ( - 180^ - 12Cli(0)2 

- lO8Cli(0 - Tvf + Cli(0)(24 - 72Cli(0 - tt)) + 72Cli(0 - tt) - 6)Ls^°^(0) 

+ ( - 5402 ^ 2167r0 - 36Cli(0)2 - 324Cli(0 - rrf + Cli(0)(72 - 216Cli(0 - vr)) 

+ 216Cli(0 - tt) - 2167r2 - 18)Ls^°^(0 - tt) + (-12Cli(0) - 36Cli(0 - tt) 

+ 12)Ls^°^(0) + (-lO8Cli(0) - 324Cli(0 - tt) + lO8)Ls^°^(0 - tt) 

+ lO80Ls^^^(0) + (2160 - 4327r)Ls^^^(0 - tt) - 4Lsf ^(0) - lOSLsf ^(0 - tt) 

- 126Lsf ^(0) - 216Lsf ^(0 - tt) + (-72Cli(0) - 216Cli(0 - tt) + 72)Lsc2,2(0) 

- lO8Lsc2,3(0) - 36Lsc3,2(0) - 3 (20^ - 126C(3)0 + 1987rC(3) + Stt^) ) (6.15) 
6.3 Four propagator integrals 

Fi\e) = \ (6.16) 
F^\Q) = -\ (6.17) 
^0(0) = Jcsc^ (^-6 + 02 + 202cii(0) + 6cos(0)-02cos(0)-4Ls^^^(0) + 40sin(0)) (6.18) 

^3^(0) = -L csc^ ( ^ ) ( - 1440 + l^-K^ - 6002 _ i2O02cii(0)2 + 144Ocos(0) + 1800^ cos(0) 
+ 1207r2 log2(2) - 1201og^(2) + 240Ls^°V)^ - 7207rLs^°V - tt) + (144OCli(0) 
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- 240(3 COS (^) - 4))Ls^'^(^) + (288OCli(0) - 1440(cos(6') - l))Ls^ ^(^ - tt) + 180Ls^^(2^) 

- 720LsJV - tt) + 144OLsLsci,i,2(0) - 2880Li4 + Ls^°^(^)(-48O0Cli(^) 

+ 1440Ls^°^(e - vr) + 24O(cos(0)6' - 9 - 2sin(0))) + 96O6'sin(0) + UmCh{9 - vr) sm(^) 
+ Lsf\e - 7r)(1440(27r - ^)Cli(^) - 720(- cos(6')e + 6 + 27rcos(^) + 2sin(0) - 27r)) 
+ 2520C(3) - 2520cos(^)C(3) - 2520 log (2)C(3) 

+ 120Cli(e) {cos{e)e^ - 26^ + Asm{e)e + 42C(3)) ) (6.19) 



F^\e) = \ (6.20) 
F-\e) = -0cot(0+^ (6.21) 

F^ie) = ^ csc2 (20 -e^-20 cos{9) + cos{e) + 8Ls^^^ {6) - 4LsP {e){e - sin(^)) 

- 120sin(^) - WCXiiO) sin(^) - 9,eCh{e - tt) sin(e) + 8Ls^°V - ^r) sin(^)) (6.22) 

Fl{e) = — csc^ ( - ) (768 - 36^2 _Qi_ 768cos(0) + 36^2^03^^) _ 96Ls^,°^(0)2 + (-192Cli(^) 
96 V^/ 

- 48(cos(^) + l))Ls^^^(^) + 96(cos(0) - l)Ls^^H^ - tt) - 96Lsi^^(^) - 48LsJ°V)(^ - sin(^)) 

- 96Lsc2,2(^)(^ - siii(^)) - 87r^ sin(0) - 480^ sin(0) + AO^ sin(0) - ASeC\i{ef sin(^) 

- 2880Cli(^ - tt) sm{e) - 960Cli(^ - vr)^ sin(^) + 96Ls_f ^ - tt) sm(0) 

+ Ls^°^(0)(-96Ls^°^(^ - tt) + 96Cli(^ - tt) sin(e) + mC\i{e){e + sin(e)) + 4S,{e + 3sin(^))) 
+ Ls^°^(e - 7r)(96Cli(0) sin(0) + 192Cli(0 - tt) sin(0) + 48(- cos(0)^ + 9 + 27r cos(0) 
+ 6sin(0) - 27r)) + Cli(0) (24 (cos(0)02 - O'^ - 6sin(0)0) - 960Cli(0 - 7r)sm(e)) - 168C(3) 
+ 168cos(6')C(3)) (6.23) 



= (6-24) 

F^^{e) = ^csc^ ( ^ ) (20cos(e) - 5cos(2^) + 16^sin(^) - 4esin(2^) - 15) (6.25) 

F^{e) = ^ csc^ ( - ) ( - 109 + 6^2 + I40cos(^) - 80^008(0) - 31cos(2^) + 29"^ cos{29) 

64 V 2 y V 
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- 32Ls^^(0) + 72^sin(^) - IGLs^^ - 7r)(4sin(0) - sin(26')) - 20^sin(2^) 
+ 8Ls^°^(^)(2e - 4sm{e) + sm{29)) + 8Cli{e){49 sm{9) - 9sm{29)) 
+ 16Cli(0 - Tr){49sm{9) -9sm{29))) 



(6.26) 



= ^ 



+ 30^2 cos(2e) + 192Ls^°^(0)2 + (384Cli(e) + 48(4cos(6') - cos(20) + 3))Ls^ ^(0) 

- 96(4cos(0) - cos(20) - 3)Ls^^V - tt) + 192Ls^^^(6') + 327r3 sin(e) + 1368esin(0) 

- 16^3 sin(^) + Lsf\9 - Tr){-A8{-Acos{9)9 + cos{29)9 + 39 + 87r cos(0) - 27rcos(2^) 
+ 18sin(0) - 5sin(2^) - 67r) - 96Cli(0)(4sin(0) - sin(20)) - 192Cli(0 - 7r)(4sin(^) 

- sin(26l))) + Ls^°^(6l)(192Ls^°^(6l - vr) - 24(66* + 20sm(6l) - 5sm(26')) 

- 96Cli(e - 7r)(4sm(e) - sm(2^)) - 96Cli(e)(2e + 4sm(^) - sin{29))) 

- 96Lsf\9 - 7r)(4 sm{9) - sm{29)) - Stt^ sm{29) - 3729 sm{29) + 49^ sm{29) 
+ 48Lsf\9){29 - Asm{9) + sm{29)) + 96Lsc2,2{9){29 - Asm{9) + sm{29)) 

+ 48Cli(6l - 7r)(186lsin(6') - 59sm{29)) + A8Cli{9fi49sm{9) - 9sm{29)) 

+ 96Cli(0 - 7rf{Aesm{9) - esm{29)) + Cli(^) (96Cli(e - TT){49sm{9) -9sm{29)) 

-24 {Acos{9)9^ - cos{29)9^ - 39^ - 20sm{9)9 + 5sm{29)9)) 

+ 504C(3) - 672 cos(0)C(3) + 168 cos(20)C(3)) (6.27) 



2880 \2j \ ' ' ' ' 

+ 2160Ls^°^(^)2 + (1440^ - 720(2^ - 47r)Cli(^) - 28807r)Ls^°\^ - tt) 

- 72OLs^°^(0 - Tvf + Lsf\9){-2880C\i{9)9 - 14400 + 288OLs^°^(0 - tt)) + 7209Lsf\9) 

- 720(0 + 7r)Lsf\9 - tt) + (432OCli(0) + 4320)Ls^^V) + (288OCli(0) - 288O)Ls^^\0 - tt) 
+ 144OLsi^^(0) + 18OLsi^^(20) - 72OLsi^^(0 - tt) + 144O0Lsc2,2(0) + 144OLsLsci,i,2(0) 




-23- 



- 2880Li4 - 720Cli(e) {9^ - 7C(3)) - 5040C(3) - 25201og(2)C(3)) (6.30) 

6.4 Five propagator integrals 

F^iQ) = -L csc^ j ( - - m'^Ch{9f + S2eCh{9)Lsf\e) - 16Lsf (^)2) (6.31) 

F^{e) = ^ csc^ - inO'^ + Ae^Ch{e) - 12Ls^^)(0)) (6.32) 

Fi{9) = -j— csc^ ( ^ ) (5537r^ + UOir^e^ + SAOtt^ log2(2) - 8401og^(2) + UA0Lsf\9f 
+ (1440 (^ - ZTT^) + 8640(27r - e)Cli(e))Ls^°V - ^r) 

+ Ls^°^(^)(lOO8OLsJ,°^(0 - tt) - 288O0Cli(^)) - 50407rLsf ^(^ - tt) + 4320Cli(e)Ls^^^(^) 
+ (17280Cli(e) - 1440z7r)Ls^^^(^ - tt) - 2880LsJ^(e) + 1260LsJ^(2^) - 5040Lsi^V - tt) 

+ lOO8OLsLsci,i,2(0) - 20160Li4 Q j - 36OzCli(0) (tt^^ + 84<(3)) 

- 2520i7rC(3) - 17640 log(2)C(3)) (6.33) 

6.5 Six propagator integrals 

F^{e) = - ^ csc^ sec (2Lsf \9)9^ - 22Ls^^\9)e - 4OLs^^^(0 - Tr)9 - 7OC(3)0 + 27Lsf ^(0) 
+ 4 {29^ - WTr9 + Qtt^) Ls^°^ {9 - tt) + 727rLs^^^ (^ - tt) + 36Lsf ^ {9 - tt) + 727rC(3)) (6.34) 

Ffo'(^) = 1 csc^ Q ((2^ - An)L4\9 - vr) - 4Ls(^V - ^) - 7C(3)) (6.35) 
F^q{9) = csc^ ^0 - 1587r^ - 2407r3^ - 75^^ - 24071^ log2(2) + 2401og^(2) 

+ (-14400 - 1440(47r - 29)Cli{9) - 144OLs^°^(0) + 288O7r)Ls^°^(0 - tt) 

+ 2880Lsf ^(0 - Trf + 1440(20 + 7r)Ls^°\0 - tt) + (2880 - 576OCli(0))Ls^^^(0 - tt) 

+ 1440LsJV) - 36OLsi^^(20) + 1440LsJV - tt) + 144O0Lsc2,2(0) 

- 288OLsLsci,i,2(0) + 5760Li4 + 5040C(3) - lOO8OCli(0)C(3) + 5040 log (2)C(3)) 

(6.36) 
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7. Virtual corrections to gg —>■ h 

In this section we will present the results for the two loop corrections to the Higgs boson production process 
via gluon fusion, with either quarks or scalar quarks running in the loops, in the region below threshold. 






Figure 4: Contributions to gg ^ h at the lowest order 
At order a^, the unrenormalized amplitude for the process gg H, is given by 

x = £EK'+>fVo(«J). 



i=sj 



4tt 



(7.1) 



where s and / denote the contributions of scalars and fermions in the loop. The Born amplitudes are given 
by the diagrams in Figure 4. They are given by 

,2 



M 



(0) 



^ab Kab Aj 



Explicit expressions for the form factors M^^^ are given in sections 7.1 and 7.2 for fermions and scalars 
respectively. Indices a and b denote the colors of the gluons, ei(pj) the corresponding polarization vectors 
and rrii the mass of the particle running in the loop. The couplings of fermions and scalars to the Higgs 
boson have been written as gfjH = ^f''^f ^i^d gssH = respectively. The constants Aj have inverse 

mass dimensions, in the case of the SM, Af = 1/v, where v is the VEV of the Higgs Boson. Finally, the 
helicity projector Kab is given by 



(0) 



(7.2) 



Kab = ea(pi) • efe(P2)pi ■P2-Pi- eb{p2)P2 ■ ea{pi 



The two loop contributions to the amplitudes can be written as 



■a. 



(7.3) 



(7.4) 
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where the infrared and ultraviolet poles have been extracted into and M-fl^-^ respectively. The form 

factors M^^gjj are finite in the limit e — 0, their explicit expressions in the region below threshold are 
given in sections 7.1 and 7.2. We have omitted contributions from diagrams with gluon self-energies in 
external lines. These drop out if one renormalizes in a heavy quark and squark decoupling scheme; the 
running of the renormalized strong coupling is then determined from the light flavors below the decoupling 
scale [40,41]. 

The singular contributions can be written in terms of the Born amplitudes. For the infrared contribu- 
tions we have 



47r V47r/x2 



+ 0{e) . (7.5) 



The ultraviolet pieces, in turn, are given by 

M^l = - (^) (2 SZ, + ^Z^,,..on ^ {Mf^)) , (7.6) 
where 6Zg is the strong coupling renormalization constant at one loop in the MS scheme, given by 

5Z, = -^(47r)^e-^-^, (7.7) 

and 52^mi, gluon are the gluonic contributions to the mass renormalization to fermions and scalar s, also in 
the MS scheme: 

-^^m^gluon = (47r)^ e"^^^ Cf - . (7.8) 

47r e 

As discussed in the introduction, the fermionic contributions have been computed in [3] and expressed 
in terms of 8 one-dimensional integrals. In turn, these integrals were computed in terms of harmonic 
polylogarithms in [25], giving the analytical result for the two loop corrections with fermions in the loop. 
Our results for these pieces fully agree with the ones quoted in [25]. 

In the case of scalars mediating the giuon-Higgs boson interaction, there are additional contributions 
originated in quartic couplings between the scalars. In supersymmetric theories, these interactions have a 
component proportional to g^. Contributions containing the quartic interactions involve, additionally, the 
mixing between different scalar quarks. Therefore, the NLO corrections associated to them contain more 
than one massive particle running in the loops. However, as the gluon couplings to the scalars are diagonal, 
the corrections due to mixing only give contributions in the form of products of one loop integrals. In what 
follows we will only consider the gluonic corrections and postpone the treatment of contributions due to 
self interactions of the scalars to a forthcoming publication. 

If the mass of the Higgs particle is significantly smaller than the mass of the particles circulating in the 
loops, the amplitudes can be approximated by their limit when — 00. These results have been obtained 
in the context of effective field theories, both for fermions and scalar loops in references [4, 42] and [27] 
respectively. Choosing the MS scheme -notice that in the infinite mass limit, the scheme dependence due 
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to mass renormalization cancels out- we obtain 



for fermion loops, whereas for scalar quarks, the amplitude is given by 

a.. 



M 



24 TT 



(7.9) 



(7.10) 



As mentioned above, this last result, which agrees with [28], does not include the self interactions of the 
scalars, thus it differs from the one quoted in [27]. Including a four scalar vertex with coupling given by 



X^a 9s i'^ij '^kl + '^il^jk ] ' ^"^^ modifying accordingly the mass renormalization of the scalars, we find 



A^s,oo,4 = ^ab Kab 



as 

24^ 



1 + 



as_ 25 
vr 6 



(7.11) 



in agreement with [27]. 

Notice that in the expressions quoted above, we have taken the limit e ^ 0, except in the infrared 
singular pieces, before evaluating the limit rrij oo. The infrared contributions contain a prefactor (m?)"'' 
that can be expanded only after combining with the real radiation pieces. 

The following two subsections contain the explicit results for the form factors M^^^^ and M-^^^ in the 
region below threshold, s < 4mf. As discussed above, in this region the natural variable to use is 6 defined 
by X = exp(z^), with < 9 < ir. In Appendix A we also give the results for the form factors as linear 
combinations of the master integrals introduced in the previous sections. 

7.1 Amplitudes for quarks 

At the one loop level, the contributions from fermion loops to the process gg ^ H when s < 4mj are 
given by 



(0) _ 1 



{4 -9^ -{4 + 9^) cos{9)) csc^ 



+e 



2 {2 77 -9) cot^ 



CSC 



Lsf (^-7r)-4cot2 



CSC 



+^ csc^ j {6-9^-29 sm{9) - 7C(3)(1 + cos(^)) - 6 cos{9)) 



40 cot { - ) csc^ ( - ) Cli{9 - tt) + 2 cot^ { - ) csc^ ( - ) Ls^°V - ^ 



+2 9 cot^ csc^ Ls^"^(0 - tt) - 4 csc^ Lsl''{9 - tt) 

-2csc^(^0 (27r-0-sin(0))Lsf (0-7r) + ^csc^(^Q (336 + 47r2 - 47r=* - 240^ 
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TT^r- (336 + 4 TT^e + TT^ (4 + ^ )) cos(0)- 144^ sin(e)- 336C(3)) 



+ 0{e^). (7.12) 



The form factor at two loops, in turn, is given by: 



(1) 

f,fin 



AT 



-AT 



y COS (^j cos{e) Lsf\e) + 18 COS ( ^ ) cos(e) Ls'l'ie - tt) 



+e COS {9 cos{e) - sm{9)) Lsf\9) 

+4 cos ((97r -5 9) cos{9) + 2 sin(0)) Ls^^)(0 - tt) 

+2 cos ((Qtt^ -10Tr9 + 29^) cos{9) + 2 {2 tt - 9) sm{9)) Lsf\9 - tt) 

- sin(^^^ ((6 + 4^2) cos(0) + 3(-2 + ^2-0 sin(^)))Cli(^) 

-1 (ll^ cos + 11^ cos (f ) + 9 sin - 7 sin ) ) LsW(^) 



+-I cosi-) (-303 + 0(4-14OC(3)) + 1447rC(3)) 



-4 sin j (18 -99^- 28 ({3) - 2 cos{9) {9 + 9^ + 14C(3))) 
- cos (^y j {9^ - 144 TT C(3) + 4 (1 + 35 C(3))) j 



CSC ,- 



16 cos^ ( - ) sin I - ] LsLsci,i,2(^) 



+y COS cos(e) Lsf ^(^) + 18 COS cos(^) Lsf\9 - tt) 

+87r cos^ (^) sin Lsf (^ - ;r) + 8 cos^ (0 sin Ls«(^) 

-2 cos^ (0 sin (0 Lsf)(2^) + 8 cos^ (0 sin LsW(^ - 

- (^16 cos^ (^^^ sin Ls^°V - tt) - ^ cos (^^^ {9 cos(e) - sin(^))^ Ls^°^(^) 

-4(8C1,(^) cos^(^) sin(^) 

-cosT^J ((97r-50) cos(0) + 2 sin(e))J Ls^^)(^-7r) 
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2 I 8 (27r - Ch{e) cos^ I ^ ) sin ( - 



cos Q j ((9 TT^ - 10 TT + 2 ^.Qg(^) + 2 (2 vr - ^) sin(0)) ^ Ls^"^ (0 - tt) 



-1 (ll^ cos (^) + 11^ COS ) - 5 sin (^) - 13 sin (|) ) Ls«(^) 
-sin(9 (-6 + 8.^-3.sin(.) + 28C(3)+cos(.)(6 + 7.^ + 28C(3))) C,(.) 

cos (0 (90^ - 4^ (1 - 35C(3)) - 1447rC(3)) 
~8 [yJ (3^' + 4^(1 + 35C(3)) - 1447rC(3)) 



-- (45 + 2 log(2)4 + 48Li4(l/2) + 42 C(3) + 42 log(2) C(3)) 
7.2 Amplitudes for scalar quarks 

The one loop form factor for scalar quarks, s < Amj is given by 
Mf)= -^(2-^2-2cos(^))csc^(Q 



csc^(-)+0(6)(7.13) 



e 

+ 7T 



4(2 TT-9) CSC* { I J Ls^°^ (0 - tt) + 8 csc^ ( | ) Ls^'^ {9 - tt) 



+ csc^ { I ) (-6 + ^2 + 6 cos(0) +20 sin(0) + 14 C(3)) 



+ 



96 



96 6* cot (^Q csc^ Ch{e - tt) - 48 csc^ Lsf\9 - Tvf 
-48 9 csc^ Lsf (0 - tt) + 96 csc^ Ls^^^ - tt) 

-48 csc^ (^^^ {-2Tr + 9 + sm{9)) Ls^^H^ - ^r) - csc^ (168 + 27r2 - 47r3 



-TT^ r - 2 (84 + TT^) cos(0) - 729 sm{9) - 168 C(3)) 



12 9' 



(7.14) 



The form factor at two loops, in turn is given by: 
r(l) _ 1 



s.fin 



iV 



— cos(0) Lsf ^ (0) + 9 cos(0) Lsf^ {9 - tt) 
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-\ sin ((-3 + Se') cos + 3 (cos (f ) - 2 ^ (^) )) Cli(0) 

~ (11 e cos(0) - 3 sin(e)) Ls^^^ (0) + ^ ^ (^ cos(e) - sin(0)) Ls5°^(e) 

+2 ((9 TT - 5 ^) cos(^) + 2 sin(e)) Ls^^^ {6 - tt) 

+((9 TT^ - 10 TT + 2 (,Qg(^) + 2 (2 TT - ^) sin(0)) Ls^°^ (0 - tt) 

(-11 ^ + - 3 cos(2 - 28 sin(0) + 17^^ gj^^^) + ^4 sij^(2 



+112 sin(^) C(3) + cos(0) (O"^ + ^ (14 - 280 C(3)) + 2887r C(3))) 



csc5 (I) 
2cos(f) 



-N 8 cos ( - J sin ) LsLsci,i,2(^) 

-— cos{e) Lsf\9) - 9 cos(0) Lsf ^(^ - vr) 
6' 



-4n cos (^0 sin (0 Lsf (^ - tt) - 4 cos sin Ls«(^) 
+ COS sin Ls^)(2 0) - 4 cos sin (0 Lsf ^(0 - n) 
~ 9 (6 cos{9) - sin(e)) LsP{9) + y (^ cos(^) - sin(e)) Ls^^^(0) 
+2 ^4Cli(0) sin(0) - ((97r -59) cos{9) + 2 sin(e))j Ls^^^(0 - tt) 

+ ( 4 (2 TT - 0) Cli (9) sm{9) + 4 sin(0) Ls^°) [9) 



-((9 TT^ - 10 TT + 2 j,Qg(^) + 2 (2 TT - 0) sin(0))^ Ls^°^ {6 - tt) 

+i sin (0 (3 cos (f ) - 6^ sin (0 - cos (0 (3 - 24^^ _ 112^(3))) Ch{9) 

1 / 1 CO 

+— I 119 + 39^ + 39 cos(2 0) + 52 sin(e) + tt^ sin(6') 

16 \ 45 

-21 9^ sm{e) + sm{e) + — vr^ log2(2) sin(0) - — log^(2) sin(0) 

3 3 

-128 Li4 ( M sin(0) - 26 sin(2^) - 112 sin(^) C(3) - 112 log(2) sin(0) C(3) 



+ cos(0) (3 e^-Ue{l- 20 C(3)) - 288 tt C(3)) 



^^^+0{e). (7.15) 
2 cos (2) 
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8. Summary 



In this paper we have computed the two-loop master integrals which are needed for the evaluation of the 
two-loop QCD amplitude in the gluon fusion process gg — > h. This is a loop induced process and generally 
requires a new calculation when modifying the particle content in the loops. 

We have automatized the evaluation of the two-loop amplitude using modern reduction methods and 
providing analytic expressions for the master integrals. We computed the master integrals using the 
differential equation method. Our results agree with the literature when a comparison is available and 
with a direct numerical evaluation which is performed with an independent method. 

In this paper we have evaluated analytically the two-loop amplitudes for gg ^ h via a quark and a 
scalar quark. The first result agrees with the result of Spira et al., in the analytic form written by Harlander 
and Kant. The amplitude for the scalar quark is a new result, and agrees with the result derived within 
the heavy squark approximation. 

The master integrals we have presented here, are relevant for other 2 — > 1 processes in the Standard 
Model and its extensions and more complicated 2^2 processes with heavy particles in the loops. 
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A. Amplitudes in terms of master integrals 

We present the results for the amplitudes in eq. (7.1) in terms of the master integrals in Section 3. We 
write the amplitudes as 



(A.l) 



A.l Fermionic amplitudes 



For the fermionic contribution at one loop, we have 




(0) 



8e{l + e)x 
(l-x)2 




(A.2) 
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At two loops: 



M 



N 



(1) ^ 



24 X (1 + x)^ - e (1 + 56 X + 238 + 56 + x''^ 



+ (9 - 360 X - 994 x^ - 360 x^ + 9 x^) + 2 (19 - 900 x - 2014 x^ - 900 x^ + 19 x^) 

+ o ^ 4 (1 - 12 X - 25 x^ + 8 x^ - 25 x^ - 12 x^ + x*^) 

x^ (1 + x)^ [ 

- 8 e (1 + 10 X + 51 x^ + 36x^ + 51 x^ + lOx^ + x*^) 



8e^ (5 + 37x + 203x^ + 214x^ + 203x'' + 37 x^ + 5x^^ 



+ 



+ 



S(l-X)2 

Cf 



s(l-x)2(l + x) 



24 (1 + xf + 20 (17 + 46 X + 17 x^) + 4 e (23 + 42 x + 23 x^) 

8 (1 - x)2 (1 - 6 X + x^) + 8 e (9 - 8 X + 30 x^ - 8 x^ + 9 x^) 



+ 8 e'-^ (21 + 8 X + 102 x^ + 8 x^ + 21 x^^ 



N 



(1-X)2 

f N 




e2 (1 - x)4 



16 X — 16 ex — 16 X 

4 s X (1 + x)^ + 2 e s X (3 + x) (1 + 3 x) + s (1 + 2 X + 122 x^ + 2 x^ + x^) 



+ 



Cf 



e(l -x)4 



lsx(l + x)^ - 4es(l - 8x - lOx^ - 8x^ + x'^^ 



N 



(l-x)4 
N 

e{l-xy 



4es(l + x)2(l-26x + x^ 



+ 



Cf 



16es(l + x) 



(l-x)2 

24 (1 + x)2 (1 - 4x + x^) - 4e (1 - 4x + x^) (7 + 26x + 7x2 
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+ 



(43 - 46x - 442 - 46x^ + 43x^) - 4e^ (1 - 4x + x^) (201 + 550x + 201x2) 

Cf 



{l-xf 



16 (1 + x)^ (1 - 4 X + x^) - 8 e (7 - 16 X - 30 x^ - 16 x^ + 7 x"") 



16e^ (10 - 19x - TOx^ - 19x^ + lOx"^^ 




N 



e2(l-x)6 



32 x^ (1 + x)^ - 16 e x^ (5 + 14 X + 5 x^) - 8 x^ (41 + 118 x + 41 x 



+ 



e(l -x) 



16s^x^ (1 + x)2 - 32es^ x^(3 + 4x + 3x^ 




N 



e(l -x)4 

N 
{l-xf 



- 8 s X (1 + x)2 + 16 e s X (1 + x^) 
32e(l + x)2 + 4(l + 6x + x2) 




Cf 



(l-x)2 
N 

(1-X)2 



8(l + 14x + x2) + 96ex 
- 80 (1 + x)2 + 2 (1 - 18x + x^) - 16 e (1 + 4x + x^) 




+ 



(l-x)2 



8(1 -6x + x^) + 8e(l + 6x + x2) + 8e2(5 + 6x + 5x 




Cf 



(l-x)2(l + x)2 



32 x^ + 16x(l + x^) - 16ex(l + 4x + x^ 




33 



N 



(1 - xY 



8sx (1 + x)^ + 24esx (1 + x) 



+ 



Cf 



{1-xY 



I6es{l-xfx + 8sx{l + xf 




Cf 



{1-xY 



N 



8s^x{l + xf (1 + 




+ 



{1-xY 

N 



4 X (1 + x)2 - 2 e X (3 + x) (1 + 3 x) 
- 96 (1 + x)2 - 32 e (5 + 16 X + 5 x^) 




I e s (1 — x)2 

- 8 (105 + 326 X + 105 x^) - 24 (161 + 446 x + 161 x^] 
16 



s(l-x)2 



4(l + x)^ + 8e(2 + 3x + 2x^) 



+ 3e2(17 + 38x + 17x2) 

A. 2 Scalar amplitudes 

The scalar contributions are given by 




+ Oie) 



(A.3) 



2e(l + e)x 
(l-x)2 




4(1 + e + e2)sx2 

0^ 



+ 



at one loop, and 




N 
es^ X 



+ 



24x-4€2(l- 108x + x^) -4e(l -24x + x2) +4e^ (5 + 462x + 5x^) 
Cf 



X (1 + x)2 



4 X (9 - 2 X + 9 x^) + 4 e (1 + 23 X + 32 x^ + 23 x^ + x"^^ 



(A.4) 
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+ 8e^ (1 + 48x + 78x2 +48j;3 + j;^^ 



+ 



+ 



S(l-X)2 

Cf 



24x-88ex-400e''x 



s(l-x)2(l + x)2 



4 (1 - (2 - 3 X + 2 x^) - 8 e (1 + 3 X + 3 + x^^ 



166^(1 + x)2(l + 3x + x2 



+ 



+ 



(1-X)2 




e2 (1 - x)4 



4x + 4ex + 4e^x 

4 s x^ - 8 e s x^ + 4 s X (1 - 10 X + x^) 



+ 



Cf 



e(l -x)4 



8 s x^ - 24 e s x^ 




+ 



+ 



(1-X)4 



24esx (1 + x) 



24 X (1 - 4 X + x^) + 40 e X (1 - 4 X + x^) + 8 X (27 - 110 X + 27 x^) 



+ 



e(l-x)4 

+ 952e^x(l-4x + x2) 
Cf 



{i-xY 



16x(l -4x + x^) + 8ex(5- 22x + 5x^) 



+ 8e2x(21 - 86x + 21x2) 
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+ 



N 



e2 (1 - x) 



32 + 96e x^ + 400 x^ 



+ 



Cf 



e(l -x) 



16 x^ + 80es^ x^ 




+ 



e{l-x) 



O 9 

> s a; — 8 e s X 




+ 



8x - 32ex 



+ 



Cf 



20x + 32ea; 




+ 



N 



{1-xy 



8x + 24ex + 80e2x 



+ 



Cf 



{1-xy 



8x — 16ex — 32e x 




+ 



+ 



Cf 



(l-x)2(l + x)2 



2 e (1 - x)2 X + 2 (1 - x)2 X - 4x (1 + x^) 




(1-x)^ 



8 s x2 - 24 e s x^ 



+ 



(1-x)^ 



4sx^ + 4esx2 




+ 



Cf 



{1-xY 



Is^ x2(l + x^) 




+ 



N 



{i-xY 



S X 



+ 



Cf 



{i-xY 



4sx2 
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(A.5) 



at two loops. 
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